A FAMILY OF CRITERIA FOR IRRATIONALITY OF EULER'S 

CONSTANT 



MARC PREVOST 

Abstract. Following earlier results of Sondow, we propose another criterion of irrational- 
ity for Euler's constant 7. It involves similar linear combinations of logarithm numbers 
L n ,m- To prove that 7 is irrational, it suffices to prove that, for some fixed m, the distance 
of d n L n ^ m (d n is the least common multiple of the n first integers ) to the set of integers 
Z docs not converge to 0. A similar result is obtained by replacing logarithms numbers 
by rational numbers: it gives a sufficient condition involving only rational numbers. Un- 
fortunately, the chaotic behavior of d n is an obstacle to verify this sufficient condition. 
All the proofs use in a large manner the theory of Pade approximation. 



1. Introduction 

In |Sondow 2003 , the author, using Beukers' integral |Beuke rs 1979J, found a criterion 
for irrationality of Euler's constant 7. It depends on the limit of the fractional part of the 
following expression 

n fc— 1 / N 2 



n re— ± / \ a 

L n = 2j2Ys( i ) {Hn-i-H^Hn + k) 

k=l i=0 ^ ' 

n 1 

where H n is the Harmonic number H n := > — . In this paper, we establish the connexion 

k 

k=l 

In \i 

between Sondow's criteria and Pade approximant of the function . Moreover, fol- 

u — 1 

lowing the same idea, we find a family of new criteria: for each integer m < n, let us set 



J n.m 



k=o ^ ' ^ ' 

d n := LCM(1, . . . ,n) and {x} = x — \_x\ the fractional part of the real x. If, for some 
integer m, if the sequence {d n ( — l) m L n m } does not converge to when n tends to infinity, 
then 7 is irrational. 

Using the property of the error term, a more precise criterion is proved here: 
if, for some integer m, the sequence {<i 2 n(— l) m L 2 n m } is asymptotically non decreasing, 
when n tends to infinity, then 7 is irrational. 
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2. SONDOW'S CRITERION WITH PADE APPROXIMANT 

Sondow considers the double integral (so-called Beukers' integral) 

r> /' ('ft -*)»(! ~yyr ixdy . 

Jo Jo (1 -xy)]n(xy) 
Applying Taylor expansion of 1/(1 — xy) around 0, he proved the following identity 

In = ( ^ ) 7 + L n - A n = 0{2-^n^) (2.1) 

where A n = I . J H n+i . After multiplication by rf 2 n, it 
i=0 ^ ^ 



arises 



( 2n \ 

d 2n [ n J 'J = d 2n (A n - L n ) + o(l). 

Sondow's criterion: 

Since d 2n A n e Z, if the sequence of fractional part ({rf2n-^n})n does not converge to 
then 7 ^ Q. 



Sebah computed this sequence for 1 < n < 2500. Its cumulative average seems to 
converge 1/2, but the mathematical proof remains to establish. 

In the following, we will show that the sequence involved in the paper by Sondow can 
be recovered by means of Pade approximation. 

Let us consider the function (ln u)/(u — l) and its Pade Approximant [n — l/n] of degree 
(n — 1/n) at the point u — 1: 

'™ _*W«) + Ru(u) (2 . 2 ) 



u - 1 £>«(«) 

where jV n and D n are polynomials of respective degree n — 1 and n, normalized by N n (l) = 
D n (l) = 1, and R n (u) = 0(u 2n ) 

From the theory of Pade approximation, it is well known that D n is related with the 
shifted Legendre Polynomial orthogonal on the interval [0,1] with respect to the Lebesgue 
weight function. Some of these expressions are 

pm = Y.( n k ) 2 t n - k {t-if (2.3) 

k=o ^ ' 



E( l)[ n t k < 2 - 4 ) 



k=0 



D n has the following expression in terms of P* 



Dn{u) = P: (1 - uY ( ^ ) 1 • (2.5) 



Replacing P* by its expressions ()2.3J2.4|) . formula (|2.5j) becomes 



2n \ 1 / n 



- ( n E V /. 



fc=0 

n—k„.k 



k=0 



The numerator N n {u) of [n — 1/n], is related with the associated polynomial of the 
denominator: 



\ / i — i ,_n \ / 



fc=l i=0 

\ — 1 n / \ / \ 1 



\n— fc+i 



Now, what is the link between Pade Approximation and Sondow's criterion? 
The definition of 7 is, primarily, 

7 = lim(if n — Inn). 

n 

An integral representation for Euler's constant is 
Formula (|2.2j) can be rewritten as 



N n {u) D n (u) R n (u)D n (u) 



In it 1 — it In it 

and 

f l l-u»N n (u) ^ | f 1 1 - «»A»(u) ^ A 1 «"^(ti)g»(«) du 

7n In It ./n 1 - U In In u 



By linearity, the second term is expanded as 



k=0 

2n \ 1 / n x 



„ i E ( I) H ^ < 2 - 8 > 

fc=0 v 7 

1 

A n . (2.9) 



2n 
n 



The first integral can be computed as following: 



ln(w) \ n J ^ ^ \ i J Jo \n(u) 

K ' v 7 fc=l i=0 v 7 " 

-1 n k-1 / s. 2 



n 



fc=l i=0 

-l 



2n , 

1 L n . (2.12) 



From the theory of Pade approximation, in formula (|2.2|) . the remainder term i? n has 
an integral representation 



D n (it) 7 l-(l-u)t 
Thanks to formulas ()2.9I2.12|) . 7 satisfies 

A n - L n r 1 u n R n {u)D n {u) ^ 



2n \ Jq ln(-u) 
n 



Thus another expression of the remainder term I n of Sondow is 



2n\ . _ r 1 u n iU«)A>(«) , 

7 - A n + L n = - / —— du 

m{u) 







o y ln(u) l-(l-u)* 
u n {l-u) 2n t n (l-t) n 



, , , ■ dt dtt 

/ 7 ln(u) (l_(l_ u )t)n+i 

thanks to integration by parts and Rodrigues formula for orthogonal polynomials. 

Thus the approximation for Euler's constant 7 ()2.1|) is a consequence of the Pade ap- 
proximation to the function (lntt)/(l — u). 



In the same manner, Pilehrood |Pilehroo5 2004J found irrationality criteria for general- 



ized Euler's constant. He defined the following linear form in logarithms 

"l m-l / v / \ 

L K,n 2 )(a) = I 1 )( fc 2 ) C^i-* + ^™2-fc - 2 ^fe) ln(m + m + a - 1) + 

m=l A:=0 ^ ' ^ ' 

£ EC-l)*" 1 - 1 /* ( * ) / ( * ~ 1 ) + m + « - 1) 

m=m+l fc=m 

Actually, following the same idea as for Sondow's criterion, it is possible to prove that 
Pilehrood' criterion comes from Pade approximations [n 2 — l,n>i] = ~ R g 2 ~^"' > (normalized 

by i? n2 _i(l) = S ni (l) = 1 to the function ln(u)/(u — 1) at the point u — 1. The linear 
form Lf ni>n2 \(a) satisfies: 



1 i i ( . ) 

3. Statement of the results 

In order to simplify Sondow's criterion, it is convenient to choose a more simple approx- 
imation. This method leads to the following theorems. 

Theorem 1. For < m < n, let us define 

Jn,m = fu^P^u) (-L + J-) du 

Jo \ mM i-wy 

mu 



then 



T -^rijm L n m -\- J n .m (^-l) 

L n , m = E(^) ( n + M(-ir +fe ln(n-m + A;+l) (3.2) 
fc=o V / V / 

Theorem 2. TTie following are equivalent: 

(a) The fractional part of d n L n m is given by {d n (— l) m L n m } = d n (—\) m J n ^ m (*) /or some 
n or m. 

f&j TTie formula (*) holds for some m and for all sufficiently large n. 
(c) Euler's constant is a rational number. 
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Figure 1. plot of values of t n = {d n (— l) m L n>m }, for m = 



A sufficient condition which involves d n L n ^ m but not J„ )m is the following 

Theorem 3. If for some integer m, {d n (— l) m L n , m } > 0.707 n infinitely often, then 7 is 
irrational. 

Computations (see Table show that this condition is satisfied for n < 1000. Numerical 
results also suggest that for each m, {c? n (— l) m L n m } is dense in the interval (0,1) and the 
cumulative average n' 1 ^]^ =1 {c/ fc (— l) m L fcim }converges to 0.5 (see Figure 1 and 2). 

To prove 7 irrational, it just suffices to show that {c? n (— l) m L n m } does not converge to 

0. 

In section |Ul we will prove the asymptotic formula 



7 = v4 n , m -L n , m + e>(4- n ). 



(3.4) 



Actually, we will prove that the error term J„ im is a totally monotone sequence (i.e. a 
sequence of moments with respect a positive measure), converging to as 4r n . 

By substituting in L n>m , ln(n + 1 + k — m) by some suitable Pade approximants, a 
sufficient condition, involving only rational numbers is the following 

Corollary 1. Let us define, for n — m + 1 = 2 p ,p G Z, 

I n \ I n + k 



L n , m ■= V [n/n) t =i + 22 



k 



(-l) n+k [n/n] t=k/[n _ m+l) 



(3.5) 
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Figure 2. plot of values of < t n >= rr 1 YTk=\ 



where [n/n] is the Pade approximant o/ln(l + t) att = 0: 

r/1 k=0 V 7 V 7 \i=0 / 

[n/n\t = ^ 

Tr^fn\(n + k\ k 

tS k )y k j 

If for some integer m, {d 2 p(— l) m I/2P+m-i,m} does not converge to when p tends to 
infinity, then 7 is irrational. 

Another sufficient condition comes from the property of the error term in the asymptotic 
formula (|3.4|) and from the upper and lower bound of the LCM(1, . . . , n): 

Corollary 2. If for some m, {g?2p(— l) m -^2P,m} is asymptotically non decreasing, then 7 is 
irrational. 



4. TWO LEMMAS 

Lemma 1. The function - — H — is a Markov- Stieltjes function. More precisely, 

m(l — u) u 

1 l f 1 l 

+ - = / w{t) dt (4.1) 



ln(l — u) u J 1 — u t 



s 



Table 1. 



{dnLnfi} {—d n L nt i} {d n L nt 2} { — <j n £ n ,3} 



n 


.7" 


{d n L n Qj 


1 


1.38868 


2 


0.56003 


3 


0.61882 


4 


2.97160 


5 


0.44808 


6 


0.31896 


7 


0.14391 


8 


0.41138 


9 


0.09667 


10 


0.06778 


11 


0.03395 


12 


0.02378 


13 


0.01719 


14 


0.01204 


15 


0.00843 


16 


0.02637 


17 


0.01639 


18 


0.01147 


19 


0.00163 


20 


0.001147 



7'' 



1.81209 
0.58439 
0.64252 
3.31151 
0.45886 
0.32064 
0.14467 
0.41543 
0.09689 
0.06781 
0.03398 
0.02379 
0.01721 
0.01204 
0.00843 
0.02637 
0.01639 
0.01147 
0.00163 
0.00114 



0.56609 
0.63428 
3.23310 
0.45719 
0.32044 
0.14460 
0.41511 
0.09687 
0.06781 
0.03398 
0.02379 
0.01720 
0.01204 
0.00843 
0.02637 
0.01639 
0.01147 
0.00163 
0.00114 



0.67030 
0.38225 
0.45913 
0.32061 
0.14465 
0.41528 
0.09688 
0.06781 
0.03398 
0.02379 
0.01720 
0.01204 
0.00843 
0.02637 
0.01639 
0.01147 
0.00163 
0.00114 



where the weight function w is 



w{t) :-- 



t (ln^(l/t-l)+7r 2 ) 



Proof. After a change of variable (u — > (1 — u) and x = 1/t — l), formula (J4.1)) is equivalent 
to 



1 1 

+ 



, ^ dt (4.2) 

m{u) l — u J x + uln x + ir 2 

The weight function w can be found with the Stieltjes inversion formula (see [Widder T 941J). 
Another way to prove formula ()4.2|) is to apply residue theorem to the function 

ft \ 1 1 
f(x) := ——- — — . 

X + U m X + 27T 

Taking the determination of In x on the complex plane cut along the positive real axis, the 
poles of / are x = —u and x = —1. 



9 



Let us define 7 r a small semi-circle z = re l9 ,—ii/2 < 9 < n/2,r > 0. the line 
z = x + ir, x running from to R, Tr the circle z = Re l9 ,0 < 9 < 2n and D~ the line 
z = x — ir, for x from R to 0. 

Now, we compute f c f(x) dx where C is the union of , Tr, D~ and j r , with the 
theorem of residue to obtain 

dx — f I — » — ■ I dx 



x + u \ In x + in Inx — inJ J x + u \ln 2 x + it 2 

, 1 1 

-227T 



lnii 1 — u 



□ 



Now, we are in position to prove a new formula for the Euler's constant 7. 
Theorem 4. The Euler's constant 7 satisfies 



ln(l +e~ z )e z 



Proof. In the integral representation of 7 (|2.fij) . let us substitute the integrand by the 
expression (|4.1j) . This leads to 

- 1 ln(l - 1) 

1 t(ln 2 (l/t- 1) 
ln(l + e~ z )e z 



7 = / g dt 



Z 2 + 7T 2 

with the change of variable t = (1 + e 2 ) -1 . 



dz 



□ 



Lemma 2. For each fixed integer m, the sequence {{~l) m J n ,m)n defined in Theorem\^is 
totally monotonic. More precisely 

(-l) m Jn,m= [ 1/ \ n Pm (v)dv 
JO 

where the weight function is 

°- I u — U , — f x 

2 



PmO) = / 7 7 — du. 

(it — It — I 7T + m 



M — U -\- V 



f 1 ( 1 1 \ 

Proof. J nm = / w n m P*(u) 1 dw appears as Legendre modified 

Jo V lriM 



moments 



of the weight function I 1 I du. 

In u 1 — u 
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For some particular cases of weight function, a sequence of modified moments can be itself 
a sequence of moments, with respect to a positive measure (see jPrevost 1994j ). Using 
Rodrigues formula for orthogonal polynomials, Lemma ^ Fubini's theorem and after n 
integrations by parts, it arises 

r-l d n / 1 1 \ 

Jn,m = / U n - mK —L^{u n {l-uT) + Z tin 

J n\ du n \mu 1 — uJ 

1 (-l) n d n f 1 1 
u n - m K {u n {l - u) n ) du / — w(t) dt 



n\ du n Jo 1 - (1 - u) t 

1 rl t l)n jn / n-m \ 



o ./o 



n! dw n V 1 - (1 - u) t 



d n 

The computation of - — I — I needs the partial decomposition of the rational 



du n \l-(l-u)t 

function — = q(u) + I — , where q is polynomial of degree 

1 — (1 — u) t \ t J 1 — (1 — u)t 

n — m — 1. 

Another expression of J n ^ m is then 

1 rl /+ _ ^\ n-m 



= I I .-(1 - »)" J (1 _ (1 _. u)f) „ +1 »W * d». 

We do the following change of variable 

u(l — u)(l — t) . , , . , u 2 — u + v 

» = i -d -,)/ e i°' ^ ~ 4 = «"> = („- „)(„-!) e ^ 

Let 4>\{v) and (/^(v) denote the two roots of the quadratic equation v = u — u 2 , 
0iO) = g ' = 2 ' 

= / V n dv (-1)- _^L_ " 1 tl — du 



'o ' J<t>i(v) \<P(v)-lJ (u- l)(u-v) 2 1 - (1 -u)(j)(v) 

which proves the lemma. 



□ 



5. Proof of Theorem [T] 

We first prove the identity (|3.1|) linking Euler's constant 7, the linear combination of 
logarithms numbers L n ^ m , the rational numbers v4 n m and the integrals J n>TO . From for- 
mula (12.61) . one substitute the integrand ( 1 J by an approximation involving 

\ln-u 1 — uj 



Legendre Polynomials as follows: 



1 (\ - u n - m P*{u) 1 - u n - m p*( u y 

In u 1 — u 

1 u n - m p:(u) + da. 











In u 1 — u 



The expression (12. 4|) of P* leads to analogous expressions L n , m . 
By linearity 

n \ n + k \ . ,. n+k / 1-u 



k=0 



^ - -e r ■ M 



7 V ^ / 7 V lnw 



= ±( n k )( n t k )(-ir +k Hn-m + k+ l) 

A„ im is treated quite differently: 

P^ satisfies the following orthogonality relation 

/ P*(u)q(u) du = 0, for all polynomial q of degree less than n. 
Jo 

I — u n-m 

Thus, by taking q(u) = , another expression for A nm is 

1 — u 

Anm= r i -™ du= r F -"»- f >'d, 



o 1 - u J 1-u 

and so A n>m is independent of < m < n. 
Let us now compute the integral in (|5.1|) . 

Legendre polynomials satisfy a three term recurrence relation which is 

(n + l)P* n+1 {u) = (2n + l)(2u - l)P n » - nP n *_» 
P * (u) = 1 P*(u)=2u-1 
Thus, v4„ jm 's also satisfy a similar recurrence relation 

(n + l)A»+l,m = ( 2 " + l)(2it - - nA»-i,n»(w) 

A = Ai = 2 

With (|5.2|) and (|5.3|) . it is not difficult to prove that 

A n ^ m = 2H n , < m < n 
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6. Proof of Theorem [2] and Theorem El 
All the arguments are based on the formula (jH.lj) . 

7 A n m L n m -\- Jfi^m ^ '* ^n7 d n A n m d n L nm -\- d n J nm 

Thus, 

dnl G Z d n L nrn — d n J n , m G Z (6-1) 

K(-l) m L n , m } = {rfn(-l) m Jn, m } (6.2) 

since A nj7n = 2H n and (—l) m J n ,m is positive. 

On the other hand, Lemma |2] implies that the sequence (J n , m ) n converges to as 4~ n . 
The numbers d n converges to infinity as e n . Thus d n (— l) m J n , m is decreasing to 0. So, for 
all sufficiently large n, {d n (-l) m J„ im } = d n (-l) m J n>m . 

If 7 is a rational number, then djv7 G Z for some A and for all n > N, d n 7 e Z. 
(c) =>- (6) follows from the previous arguments. Thus (c) =>- (6) =>■ (a) => (c). 

In f!6.2|) . we substitute d n by an upper bound : d n < e im9n |Rosser et al. 1962] . 
Thus Vn, d n {-l) m J^ m < e L039 n 4~ n < 0.707" and Theorem |3] is proved. 

7. Proof of Corollaries 
1) In the numerical computation of formula 



J n,m 

k=0 



E(j)("£*)(-ir*M»-m + * + i>, 



the problem is the evaluation of logarithmic functions. 

A mean to avoid this drawback is the substitution of ln(n — m + k + 1) by some suitable 
approximations, enough good to keep the irrationality criteria. We will show now that 
Pade approximants satisfy this condition: 
another expression of L nm is 



L n>m = ln(n -m + l) + jr( n k )( n+ k k ) ^ n+ " ln ( X + 



k 



k=0 

The Pade error for the logarithmic function is 



n — m + 1 



(_i\n n+l pi pn /i _ f)n 



Let us set 



L' n , m := ln(n-m+l)+]C(2) ("**) 



(-l) n+fe [n/n] i=fc/(n „ m+ i) 



We have to evaluate the difference b n ^ m := L n ^ m — L' n>m . For sake of simplicity, we set 
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E 

fc=0 
n 

E 

k=0 



n + k 
k 

n + k 
k 



(-l)^(ln(l + G)-[nM =c J 



r 

T) s/l ' 



n+l 



1 j.n 



r (l -t) n 



Since Cfc G [0, 1] and P* has all its roots in [0, 1], 
hand, the integral 



An 
Sfc 



1 -in 



t n (1 -t) n 



(1 + Cfc*) 



n+l 



< 4" 



< 4" 



1 



< 



. On the other 



i 



W Ck 



So, 



|<$n,m| < ^ 



k=0 
n 

£ E 

fc=0 

< 



n + k 
k 

n + k 
k 



n+l 



l^n*(-l)| 4 ^Ck 



1 t n (1 - t) n 
o (1 + Ck 

1 



df 



|nP*(-l)| 



4-E 



fc=0 



n + k 
k 



\nP*(-l)\ 



4- n ip:(-i)i = (^4 



n\-l 



The goal is partly reached since the error between L n>rn and its approximation is less 
than J nm . Now, let us consider the approximation of ln(n — m + 1). It is difficult to 
approximate this number (which tends to infinity) with an error lees than 4~ n . So, we 
consider sequences of integers n, such that n — m + 1 is a power of 2: n — m + 1 = TP . 
With this hypothesis, ln(n — m + 1) = p In 2 

(—1)" [ l t n (l — t) n 

In (|7.1j) . if x = 1, In 2 — [n/n] x=1 = — — / — — ; — t^jT dt. The asymptotics for 



P*(-l)io (1+ 
Legendre polynomials are well known 

P„(a) ~ (a + v/a^T)™, for a G R \ [-1, 1]. 
Thus shifted Legendre Polynomials satisfy 



P*(t) ~ ((2* - 1) + 2 v / t^ r t) n , for t G R \ [0, 1]. 
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The maximum of the fraction f ^ — ^ ^ for t G [0, 1] is obtained for t = y/2 — 1, and 
its value is (3 — 2^/2). Thus 

|ln2- [n/nL-A < (3 - 2v f )n In 2 

For n - m + 1 = 2 P , ln(2 p ) - p [n/nj^i < p (3 - 2^2f n which is a o(4~ n /ra). At last, the 



error 



satisfies 



< {4T n /n) 



and the corollary 1 is proved. 

2) For the proof of Corollary |2l we exploit the property of totally monotonic sequences 
(TMS). 

A sequence u n is called TMS if there exists a non negative measure d[i with infinitely many 
points of increase such that 



j 1 oo 

Wn G N, u n = / x n d(j,(x). 
Jo 



If the support of the measure d\i is the interval [0, l/R], then Wn,u n+ i/u n < R and 
lim n+ = R. If R — 1, it is equivalent to 

n U n 

\/n G N, VA; G N, (-l) k A k (u n ) > 

where A°(u n ) := u n and A k+1 u n = A k u n+1 - A k u n . (see jWidder 194lj . p. 108). 

The previous properties can be applied to the sequence J n>m for which we prove some 
convergence properties. If they are not satisfied by {d n (— l) m L n , m } then 7 is irrational. 

First we will prove that J„ im satisfies dm{—^) m Jin,m < d n (— l) m J„ im : the numbers d n 
and J„ im satisfy 2 n < d n < e 1 039 n (see |Tenenbaum 1990] , p. 12-13 for the lower bound 



and [Rosser et al. 1962] for the upper one) " +1 ' m < 1/4 (property of totally monotonic 

J ri 

sequence |Widder 1941j . p. 135). 

dnJn, 



'n,m 



> ; ^ = 4" > 1.0014 



dlnJln,m & 

Thus, for all integer m, (c?2p( — l) m ^2P,m) peN is a positive decreasing sequence, converging 
to 0. So, if ({rf 2 p(— l) m -^2P,m}) p is non decreasing for p greater than any integer, then 7 is 
irrational. 



I would thank my colleague S. Eliahou for reference ( |Pilehrood 2 004J) 
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